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Abstract 

We consider deeply virtual Compton scattering on a photon target, in the generalized Bjorken 
limit, at the Born order and in the leading logarithmic approximation. We interpret the result 
as a factorized amplitude of a hard process described by handbag diagrams and anomalous gen- 
eralized parton distributions in the photon. This anomalous part, with its characteristic In(Q^) 
dependence, is present both in the DGLAP and in the ERBL regions. As a consequence, these 
generalized parton distributions of the photon obey DGLAP-ERBL evolution equations with an 
inhomogeneous term. 



1 Introduction 

The photon is a fascinating object for QCD studies. Among its many aspects, its parton distributions 
have been the subject of much work since the seminal paper by Witten [T]. The mixing of the non 
local electromagnetic operator F'^'^{z)F'^''{Q) with the ip{z)ip{0) quark-antiquark one, on the light- 
cone, yields quite unique features in the analysis of the imaginary part of the forward 7*7 — > 7*7 
amplitude. A renormalization group analysis allows then to define a parton distribution in the photon 
which factorizes in the structure functions measured in inclusive deep inelastic scattering (DIS). The 
pointlike nature of the photon enables one to fully determine its leading expression in InQ^, through 
a Born order calculation in aem and an inhomogeneous QCD evolution equation which can be solved 
in the leading logarithmic approximation without assuming an additional initial condition. 

In this work we uncover an analogous structure in the case of the generalization of the Bjorken 
scaling regime of exclusive hard reactions such as the amplitude for deeply virtual Compton scattering 
(DVCS), 7*7 77. At first sight, this looks imrealistic since QED gauge invariance demands that 
all diagrams contributing to the DVCS amplitude be considered and not all of them have a topology 
compatible with a partonic interpretation based on the QCD collinear factorization of the scatter- 
ing amplitude. We show nevertheless how to define the anomalous generalized parton distributions 
(GPDs) in the photon. It can serve as a basis for a reliable perturbative calculation of a GPD, which 
is of utmost importance to test various features of GPDs such as sum rules, crossing properties [2] 
with generalized distribution amplitudes 3J and positivity limits [3]. 

Moreover, the parton distributions in the photon have turned out to be of experimental importance 
in a number of accessible processes, both in e^e^ annihilation and in photoproduction. In the same 
spirit and even if DVCS on the photon seems to be more a subject for an academic study than for a 
phenomenological analysis of forthcoming data, it may turn out that other exclusive reactions with 
photon GPDs (e.g. 7*7 pj) are feasible. Let us stress, nevertheless, that the phenomenological use 
of the photon GPDs requires to include non-pointlike, hadronic contributions, which effectively goes 
beyond the leading logarithmic approximation considered below. 

2 The DVCS process 

Deeply virtual Gompton scattering on a photon target 

i*{qh{pi) i{q')i{p2) (1) 

involves, at leading order in aem, and zeroth order in as the six Feynman diagrams of Fig. [1] with 
quarks in the loop. Our aim in this section is to calculate this amplitude and present it in the form of 




Figure 1: The Born order diagrams for 7* 7 ^ 7 7 

an integral over the quark momentum fraction x. Since we found this calculation rather instructive, 
we shall describe with some details how we perform it on the example of the box diagram of Fig. lA 



1 



and then exhibit the total result coming from the sum of the six diagrams. For simplicity, we restrict 
in this short paper to the case where At, the transverse part of A = p2 — Pi, vanishes, that is zero 
scattering angle (but still non- forward kinematics). We expect to generalize the kinematics in a future 
work. Our conventions for the kinematics are the following: 

q = -'2£,P + n, q' = n, 

Pi = {1+0p, P2 =Pi + A = . 

where p and n are two light-cone vectors and 2p ■ n = s = ^ . 
The momentum k in the quark loop is chosen to be 

fc^ = (a; + + /3n^ + kr , 

with k^ = — k^. 

The amplitude of the process 7*7 — > 77 with one virtual and three real photons can be written as 

A = e^e:;ei„e;^T'-"^ (2) 

where in our kinematics the four photon polarization vectors e{q), e'*{q'), ei(pi) and £2(^2) are trans- 
verse. The tensorial decomposition of T'^'^"^ reads [5] 

r'"'"''(AT - 0) ^ Ig^'^gfw^ + l [g^'^gf + ^^fff " 5^5?^) + \ (gi^.'^gf - glfg^^) W, , (3) 

and it involves three scalar functions Wi, i — 1, 2, 3. 

The integration over k is performed as usual within the Sudakov representation, using 

d-^k = ^dxdf3(fkT = '^dxdl3d\f? 

Let us present in some details our calculations for diagram A, and afterwards only give final results 
for the other diagrams. We start with the straightforward calculation of the imaginary part of Wi 
for this diagram. Performing the dj3 and dk^ integrations with the help of the (5-functions defined by 
Cutkosky rules, we get for the symmetric Wi projection g^ g^^ '■ 

where the x— integration range is restricted by the positivity conditions attached to the (5-functions, 
to 



C + — < .T < 1 . (5) 

s s 



This leads for this diagram to 



Dtsc - 2ilm = - l^Nc log ^ , (6) 



where we restricted ourselves to the logarithmic terms and traded s for in the argument of the 
logarithm. 

The corresponding result for the discontinuity of W3 attached to the ^ ^ v antisymmetric tensorial 
structure is 



Disc = 2ilm = '-e\Nc log ^ , (7) 

and the third discontinuity Disc W-^ vanishes. 

The calculation of the real part TZe is more involved. We get for TZe 

TZei^^^ [ dxd/Sdk^ TrAi 
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where TrAi ~ Tr{'^^{k + q + m)j^{k + m)"f^{k — pi + m)j^{k + q — h + m) . 

We first integrate in (3 using the Cauchy theorem. The propagators induce poles in the complex 
/9-plane with values 

+ — irj + — irj 

k^ + — irj k^ + — i?7 , 

Since the four poles are all below the real axis for a; > 1 and all above the real axis for x < — ^, the 
only region where the amplitude may not vanish is 1 > a; > — ^. In this way we identify two different 
regions : 

• the DGLAP region where ^ < x < 1, for which one may close the contour in the upper half 
plane and get the contribution of the pole (3^ , 

• the ERBL region where — < a: < for which one may close the contour in the lower half plane 
and get the contribution of the pole /3i . 

A number of technical simplifications are now helpful. The integral over k^ in Eq. ([5]) is UV divergent. 
On the other hand, it is a well-known classical result of QED that the sum of integrals corresponding 
to the six diagrams of Fig. 1 is finite, so we will separate UV divergent terms of each diagram in an 
algebraic way. One easily verifies that, at the leading logarithmic approximation we are interested in, 
the trace may be simplified by taking the limit — *■ 0. The integral in Eq. ([8]) may then be written 
as 

I- = I- + If ^ 2.. r dx f dl^^IlM^ _ 2.. f dx f dU^^^^f-^^^ , (10) 
with the pole values of the propagators : 



s[\^ + m^fjl - - + m^) + s{l - - x)) 

-2s(k^ + m^Yi{\ + g)(2e(k2 + m^) + sj^ ~ x'^)) 



The traces are simple polynomials in k^, which may be written as {i = 1,4) TrAi{P = jSi) ^ 
a^k^ + 7i(5^k^ + SiQ'^. Power counting in k^ shows that the a^k"' term in these integrals is ultraviolet 
divergent. Our aim is to recover the UV finiteness of the sum of diagrams before performing the 
integration over x. For this, we simply separate the ultraviolet divergent part of the integrand, . 
This yields for the diagram A and the DGLAP region ^ < a; < 1 (at the leading logarithm level) 

It = 2^. C dx I dk^Zl^g^ + C dx [ dk'^ , (11) 

with the regularized numerator TrA\{[3 — P4) — 74Q^k^ + 54^Q^. Simple algebra leads then to: 

The corresponding result for the ERBL region is 

If = t d. (- + f](^-^^^- + ^ og^ + 2.. r dx jdk^^. (13) 

The same procedure is applied to the five other diagrams of Fig. 1. Diagrams B and C contribute 
to the same regions ^ < a; < 1 and — ^ < a; < ^ as diagram A. Diagrams D, E, F contribute to 
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the "antiquark" regions — 1 < a; < — ^ and — ^ < a; < ^. The results for diagrams B, D and E are 

similar to the one for diagram A. It is worth emphasizing that diagrams C and F in our leading 
logarithmic approximation contribute only to the ultraviolet divergent contribution and not to any 

2 

log ^ terms. This fact will lead to a handbag dominance interpretation of the leading logarithmic 
result, as discussed later on. 

At this point, we observe that the ultraviolet divergent contributions cancel out, before integration 
in the x variable, and moreover separately in the two sets of diagrams (A, B, C) and (D, E, F). 

The final result of our calculation of the real part of the DVCS amplitude reads 



UeWi 



Nc 



2 



dx 



2x 



+ e{^-x)e{i + x 



x{l-i) 



-e{-x-o- 



+ (1 - xf - e 

(1 + xf - e 



'^(1+0 - - - 1 

The real parts of W2 and W3 are calculated in the same way. We get 

'R.eW2 = 



2 



(14) 



(15) 



and 



2 7r2 



dx 

1 + 



e{x-o 



- - x)6{^ + x) 



1-g 
1 + C 



+ e{-x-^) 



- (1 + xf - e 
i-e 



m 



In^. (16) 



The real and imaginary parts of the amplitudes may be put together in the simple form 



Wi 



etNc 



dx 



2x 



+ 



2 7r2 

- x)e{^ + x) 



{x - + iri){x + ^ — irj) 
x{l-0 



0{x-O 



a;2 + (1 - xf - e 



0{-x-O 



x^ + (\+xf 



and 



W3 



e^Nc 



Nc r 



dx 



2^ 



(a; - ^ + iri){x + ^-iri) 



e{x-i)- 



2 

, m 



{i-xf-e 



- ei^-x)ei^ + x) 



9{-x-0 



x^-{i+ xf - e 
1-e 



In 



m 



(17) 



(18) 



We now want to interpret this result from the point of view of QCD factorization based on the 
operator product expansion, yet still in the zeroth order of the QCD coupling constant and in the 
leading logarithmic approximation. For this, we write for any function J^{x, ^) the obvious identity : 



^{x, log g2 = -^V^' log jjp" + •^(^' log -g2" 



(19) 



where Mp is an arbitrary factorization scale. Wc will show below that the first term with log may 



be identified with the quark content of the photon, whereas the second term with log -gf corresponds 
to the so-called photon content of the photon, coming from the matrix element of the two photon 
correlator Ap(— |)A,y(|) which contributes at the same order in a^m as the quark correlator to the 
scattering amplitude. Choosing Mp = Q"^ will allow to express the DVCS amplitude only in terms of 
the quark content of the photon. 
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3 QCD factorization of the DVCS amplitude on the photon 



To understand the results of Eqs. (|17m8p within the QCD factorization, we first consider two quark 
non local correlators on the light cone and their matrix elements between real photon states : 



Ft = 



dz 
2^ 



e^-^{^{p')\q{--N)^.Nq{-N)\^{p)) 



and 



'-{^{j>')\q{--N)^.N^\{-N)\^{p)) 



(20) 



(21) 



where we note N — n/n.p and where we neglected, for simplicity of notation, both the electromagnetic 
and the gluonic Wilson lines. 

There also exists the photon correlator f'{~^N)F''^ {'^N) (where F^^' = iV^F'^^), which mix 
with the quark operators [T], but contrarily to the quark correlator matrix element, the photonic one 
begins at order ag,„, as seen for instance in the symmetric case : 



'^e-^ia{p^)\F^^{-^-N)F^^{^N)gTMvr)) 



-.grep(pi)<(P2)(i - e)m + ^) + -5(1 - ^)\ ■ 

(22) 

The quark correlator matrix elements, calculated in the lowest order of a^m and as, suffer from 
ultraviolet divergences, which we regulate through the usual dimensional regularization procedure, 
with d = 4 + 2e. We obtain (with 4 = i + 7^; - log47r) 



log TO'' 



with 



5 B{\ > x> i) 



1 



1 



-i>x> -1) 



(23) 



-^^(e > - > -0 



(24) 

for the [I V symmetric (polarization averaged) part. The corresponding results for the antisymmetric 
(polarized) part read 



- + log mf 



F(x.i) 



with 

F{x,i) 



(1 



1 



(1 > a; > 



(\ + xf~e 



-i>x> -1) 



(25) 



(26) 

Let us stress again that here we concentrate only on the leading logarithmic behaviour and thus focus 
on the divergent parts and their associated logarithmic functions. The ultraviolet divergent parts are 
removed through the renormalization procedure involving both quark and photon correlators (see for 
example [S]). The 4 terms in Eqs. (HSHH]) define the non-diagonal element Zqq ff oi the multiplica- 
tive matrix of renormalization constants Z. The 4 terms are then subtracted by the renormalization 
of quark operators. The renormalization procedure introduces a renormalization scale which we here 
identify with a factorization scale Mf in the factorized form of the amplitude (see Eq. (l30|) below). 
Imposing the renormalization condition that the renormalized quark correlator matrix element van- 
ishes when the factorization scale Mf = to, we get from Eq. (|23p for the renormalized matrix element 



= ^|^5rep(Pi)ei^(P2)log-^F(x,^) 

F 



and a similar result for the antisymmetric case of Eq. (^51 

Nce^ 



F'n 



(27) 



(28) 
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Eqs. (P71 [^5)1 permit us to define the generalized quark distributions in the photon, H'^{x,£,,0), as 



pi 



and to write the quark contribution to the DVCS amphtude as a convolution of coefficient functions 
and distributions Hf 



W? = I dxCl{x)Hl{x,i,Q) 



dxCl{x)Hlix,^,0) 



(30) 



where the Born order coefficient functions Cy^j^ attached to the quark-antiquark symmetric and 
antisymmetric correlators are the usual hard process amplitudes : 

1 1 



^VjA — 



± 



(31) 



We recover in that way the In term in the right hand side of Eq. (|19p. 

The photon operator contribution to the DVCS amplitude at the order considered here, 
involves a new coefficient function calculated at the factorization scale Mp, which plays the role of 
the infrared cutoff, convoluted with the photon correlator ([^^ or with its antisymmetric counterpart. 
The results of these convolutions effectively coincide with the amplitudes calculated in Section 2 with 
the quark mass replaced by the factorization scale m — > Mp, and leads to the second term in the right 
hand side of Eq. ([TO]) . The triviality of Eq. in fact hides the more general independence of the 
scattering amplitude on the choice of the scale Mp which is controlled by the renormalization group 
equation. 

We still have the freedom to fix the factorization scale Af|. in any convenient way. Choosing 
kills the logarithmic terms coming from the photon correlator, so that the DVCS amplitude 
is written (at least, in the leading logarithmic approximation) solely in terms of the quark correlator, 
recovering a partonic interpretation of the process, see e.g. 0. 

4 The GPDs of the photon and their QCD evolution equations 

We have thus demonstrated that it is legitimate to define the Born order generalized quark distribu- 
tions in the photon at zero Ay as 



Ml 



mix, to) = 



Ncel 



(1 



9ii~x)e{^ + x) 



Hl{x,tO) 



47r2 



WW) ''^ 

2 - (1 - a;)2 



-e 



+ (1 + xf 



e{^^x)e{i + x)^—^^ +6{-x-s,)- 



i-e 



{i+x)^-e 



In 



(32) 



In ^-7T 



Since we focussed on the logarithmic factors, we only obtained the anomalous part of these GPDs. 
Their x— and ^—dependence are shown on Figs. 2 and 3. The <^ — > limit of Eqs. ([5^ leads to 
the usual anomalous Born order quark content of the photon. The Hf{x,£,,0) GPDs are continuous 
functions of x at the points x = but their derivatives are not. 

On can verify that these GPDs obey the usual polynomiality conditions 



/ 



1 



1 



dxx'^P+^ Hf{x,tO) 



p + 1 

2 



(2p + 3)(p + 2) 
2 



3p- 



e2p+2 



P- 
p 



m=0 



2m 



(34) 
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the other moments vanish. 




X 



0.2 0.4 0.6 0.8 



Figure 2: The unpolarized anomalous photon GPD / {Nc e^/ {Atv'^) In^) at Born order and At = for 
^ = (full line), 0.1 (long dashed line), 0.3 (middle dashed line), 0.5 (small dashed line). 




Figure 3: The polarized anomalous photon GPD H§/(Nc eq/{'iTY^) In £5-) at Born order and At = for 
5 = (full line), 0.1 (long dashed line), 0.3 (middle dashed line) 0.5 (small dashed line). 



Switching on QCD, the photon GPDs evolve from this QED-like parton model result through 
modified DGLAP-ERBL evolution equations 



1 



(35) 



where ff{x,^) are defined by the r.h.s. ofEqs. (|32l[55)) . as the corresponding functions which multiply 
InQ^/m^. The QCD kernel Vns is well-known [5] to the leading order accuracy 



Vns{x,x') = -r-Cp 

4 TT 



p{x,x') 



l+X 

l + x' 



+ \x ^ —X, x' 



pix, x') = 9{x' >x>-l)- e{x' <x< -1) 



(36) 



with Cf = (A^c - l)/(2A^c)- As in the forward case, the presence of the inhomogeneous term allows 
to fully evaluate the GPDs. We will address this point in a future work. 
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5 Conclusions 

We derived the leading amplitude of the DVCS (polarization averaged or polarized) process on a 
photon target. We have shown that the amplitude coefficients factorize in the forms shown 
in Eq. (|30p . irrespectively of the fact that the handbag diagram interpretation appears only after 
cancellation of UV divergencies in the scattering amplitude. We have shown that the objects Hf{x, t) 
are matrix elements of non-local quark operators on the light cone, and that they have an anomalous 
component which is proportional to \og{Q^/m^). They thus have all the properties attached to 
generalized parton distributions, and they obey non-homogeneous DGLAP-ERBL evolution equations. 

The extension of the results of this work beyond the leading logarithmic approximation requires 
taking into account the non-pointlike, hadronic contribution of the photon wave function. In the 
diagonal case of the photon structure function, it is done in a model dependent way, by invoking 
the vector meson dominance hypothesis (see second Ref. 1 ). This approach encounters, however, 
well known difficulties related to the impossibility to separate in a precise way those non-pointlike 
contributions from the next to leading order corrections to the processes with the pointlike coupling of 
the photon 9J. In the non-diagonal case considered here, we may apply a similar strategy and define 
non-pointlike photon GPDs by relating them to the vector meson GPDs or to the 7 —> vector meson 
transition distribution amplitude [TUj . 
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